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Abstract 

We study the approximation of stochastic differential equations driven by a frac- 
tional Brownian motion with Hurst parameter H > 1/2. For the mean-square error 
at a single point we derive the optimal rate of convergence that can be achieved 
by any approximation method using an equidistant discretization of the driving 
fractional Brownian motion. We find that there are mainly two cases: either the 
solution can be approximated perfectly or the best possible rate of convergence is 
ra - ^ -1 / 2 , where n denotes the number of evaluations of the fractional Brownian mo- 
tion. In addition, we present an implementable approximation scheme that obtains 
the optimal rate of convergence in the latter case. 
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1 Introduction 



In this article, we study the pathwise approximation of the equation 

dX t = a(X t ) dt + <r(X t )dB t , t e [0, 1], (1) 
Xq — x 0) 
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where Xq G M and (Bt)te[o,i\ is a fractional Brownian motion with Hurst param- 
eter H G (1/2, 1). Equation (1) is understood as a pathwise Riemann-Stieltjes 
integral equation, see e.g. [27], [31] and [37]. Recent applications of stochastic 
differential equations driven by fractional Brownian motion include e.g. the 
noise simulation in electronic circuits ([11]), the modelling of the subdiffusion 
of a protein molecule ([16]) and the pricing of weather derivatives ([3], [4]). 

The type of approximation methods we are concerned with produce an ap- 
proximation to X 1 using an equidistant discretization of the driving fractional 
Brownian motion, that is 

B\/ n , B 2 / n , • • • , B\. 

The error of such an approximation method will be measured by the mean 
square norm (E| • I 2 ) 1 / 2 . Thus, the best possible approximation method of the 
above type is clearly the conditional expectation 

— E( Xi | B 1 / n , B 2 / n , ...,B ± ). 

For stochastic differential equations driven by Brownian motion, the optimal 
pointwise approximation (in the mean-square sense) of the solution is a well 
studied problem, also for non-equidistant and adaptive discretizations. See, 
e.g., [7], [24], [25], [6], [5] and [20]. However, for stochastic differential equations 
driven by fractional Brownian motion there are only a few known results for 
mean-square approximation, mainly for equations with additive noise ([21], 
[22]) or with a linear diffusion coefficient ([19]). In [21] and [22] the Euler and a 
Wagner-Platen-type method for equations with additive noise are studied and 
their exact rates of convergence are given, while in [19] the convergence order of 
an Euler-type method for a quasi-linear Skorohod-type equation is determined. 
Moreover, the asymptotic error distribution of several approximation schemes 
for equation (1) is derived in [23] and [13]. 

Throughout this article, we will impose rather strong assumptions on the 
drift- and diffusion coefficients, which has technical reasons. See Remark (i) 
in Section 5.1 for a detailed discussion. In particular a and a are supposed 
to be bounded and a is also strictly positive. Nevertheless, we think that this 
article will give a lot of structural insight in the approximation of stochastic 
differential equations driven by fractional Brownian motion. 

If the drift and diffusion coefficient commute, i.e. a' a — aa' = 0, then it is easy 
to show that X 1 does not depend on the whole process (S t ) te [ 0)1 ] but only on 
B\. This implies that 

(E|X 1 -XS n) | 2 ) 1 / 2 = 0. (2) 
More generally, we can show the following relation: 
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There exists a mapping g e C 2 ([0, 1] x R; R) such that X t = g(t, B t ) holds 
for all t E [0, 1] almost surely, if and only if (a'a — aa')(x) = for all iGl 

This extends a particular case of a well known result for stochastic differential 
equations driven by Brownian motion. See e.g. [36]. 

Otherwise, if a and a do not commute, we have the following upper and lower 
error bound for the error of the conditional expectation 

liminf n H+l / 2 (E|Xi -X[ n) | 2 ) 1/2 > a H (jf* \Ey t \ 2 dt^ (3) 

and 

limsup n H+1 / 2 (E|X! -xf^ 2 ) 1 ' 2 < (3 H ( [* E\y t \ 2 dt) , (4) 

n->oo VO / 

where ajj and (3h are two constants, depending only on H and the random 
weight function (^) t6 [o,i) is given by 

y t = {aa> - a'a)(X t ) exp Qf* a'(X s ) ds + jf* a'(X s ) dB^j , t £ [0, 1]. 

If the constant on the right hand side of equation (3) does not vanish, then 
the conditional expectation has exact rate of convergence n~ H ~ 1 ^ 2 . This is 
satisfied for example, if (acr' — a'a)(x ) ^ 0. Consequently, in this case there 
is no approximation method using an equidistant discretization of the driving 
fractional Brownian motion that can obtain a better rate of convergence than 

n -H-l/2 m 

The conditional expectation is clearly in general not an implementable method 
for the approximation of stochastic differential equations driven by fractional 
Brownian motion. Therefore, we also consider here an extension of McShane's 
method for Stratonovich SDEs, see [17] and [20] for a related scheme for Ito 
SDEs. Our extension of the McShane method is defined by X^ = x and 

*Si = *l n) + «(*J°)A + o{X^)A k B + l -aa\X^){A k Bf 
+ h aa > + a 'a)(xl n) )A k BA + \aa\X { k n) )A 2 
+ 1 -(a 2 a'' + a(a') 2 )(xi n) )(^Bf 

for k — 0, . . . , n — 1, where A = \jn and A k B = -B(fc+i)A — B k &. 
We show that 

Jim o n^ +1 / 2 (E|X 1 -)CA")| 2 ) 1 / 2 = fa E\y t \ 2 dt) ^ . (5) 
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Hence this implementable approximation method has exact rate of conver- 
gence n~ H ~ 1 / 2 '', if a and a do not commute, and thus obtains the same con- 
vergence rate as the conditional expectation in this case. 

The article is structured as follows. In the next section we recall some facts 
about fractional Brownian motion and stochastic differential equations driven 
by fractional Brownian motion. In the Sections 3 and 4, we state and prove our 
results for the error of the conditional expectation, while McShane's method 
is considered in section 5. A technical proof of an auxiliary result is postponed 
to the Appendix. 



2 Preliminaries 

2. 1 Fractional Brownian motion 

Let B = (Bt)te[o,i\ be a fractional Brownian motion with Hurst parameter 
H G (0,1) defined on a complete probability space P), i.e., B is a 

continuous centered Gaussian process with covariance function 

R H (s, t) = l -{s 2H + t 2H -\t- s\ 2H ), s, t G [0, 1]. 

For H = 1/2, B is a standard Brownian motion, while for H ^ 1/2, it is 
neither a semimartingale nor a Markov process. Moreover, it holds 

(E\B t -B a \ 2 )^ = \t-s\ H , Me [0,1], 

and almost all sample paths of B are Holder continuous of any order A G 
(0,H). 

Let us give a few facts about the Gaussian structure of fBm for H > 1/2 and 
its Malliavin derivative process, following [29] and Chapters 1.2 and 5.2 in 
[30]. Let 8 be the set of step functions on [0, 1] and consider the Hilbert space 
H defined as the closure of £ with respect to the scalar product 

(l[o,t],l[o,«])w = Rii(t,s), s,te [0,1]. 

The mapping l[ 0)t ] i— > B t can be extended to an isometry between Ti and its 
associated Gaussian space. This isometry will be denoted by (p i— > B (</?). Note 
that 

(<P,p)n = [\K* H <p)(s)(K* H p)(s)ds : <p, P e£, (6) 

J 
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where 

- 1 , ,dK 



(*»(s) = J tp(r)-^-(r, s) dr, up G £, s G [0, 1], 



with 



K H (t, s) = c H s^- H jf (« - s) H -^u H -^ du l m (s) 

and 

/ H(2H-1) \ 1/2 
° H ~ \P(2-2H,H-l/2) ) ' 

Here (3 denotes the Beta function. Moreover, we have L 1 / // ([0, l])cH and in 
particular 

{ip,ij) H = a H f f (p(r)^(u)\r-u\ 2H ~ 2 drdu (7) 
Jo Jo 

for (p,ip G L 1 / // ([0, 1]) with a H = H(2H — 1). (For a characterization of the 
space Ti in terms of distributions see [15].) 

Let / : R" — > K. be a smooth function with compact support and consider 
the random variable F = f(B tl , . . . , B tn ) with ti G [0, 1] for % — 1, . . . , n. The 
derivative process of F is the element of L 2 (fl; Ti) defined by 

D S F = ^-^(B tl ,...,B tn )l [0:ti] (s), se [0,1]. 

In particular D s B t = l[o,t](s). As usual, D 1 ' 2 is the closure of the set of smooth 
random variables with respect to the norm 

||F|| 2 2 = E|F| 2 + E||DF||^. 



If Fi, F 2 G D 1,2 such that F 1 and ||DFi|| w are bounded, then F X F 2 G D 1,2 and 

DF 1 F 2 = F 2 DF 1 + F 1 DF 2 (8) 

Moreover, recall also the following chain rule: For F G D 1,2 and g G C 1 (IR) 
with bounded derivative we have g(F) G D 1 ' 2 and 

Dg{F) = g'(F) DF. (9) 

The divergence operator 5 is the adjoint of the derivative operator. If a ran- 
dom variable u G L 2 (Q;Ti) belongs to the domain dom(5) of the divergence 
operator denoted, then 8(u) is defined by the duality relationship 

E(F5(u)) = E(DF,u) n 
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for every F G D 1 ' 2 . 



If u = (ut)te[o,i\ is a stochastic process with Holder continuous sample paths 
of order A > 1 — H, then the Riemann-Stieltjes integral with respect to B is 
well defined. If u moreover satisfies u t G IB) 1 ' 2 for all t G [0, 1] and 

sup E|w s | 2 + sup Fi\D r u s \ 2 < oo, 

se[0,l] r,se[Q,l] 

then we have u G dom(5) and the relation 

f 1 u t dB t = 5(u) + a H [ f 1 D s u t \s - t\ 2H ~ 2 dsdt (10) 
Jo Jo Jo 

holds. For the Skorohod integral of the process u we have the isometry 

E\5(u)\ 2 = a H [ Eu s u r .\s-r\ 2H - 2 drds (11) 
J[o,i] 2 

+ a 2 H ( ED r u s D r ,u s/ \r - s'\ 2H ~ 2 \r' - s\ 2H ~ 2 drdr' dsds' . 



In what follows, we will also the require the Wiener-Chaos decomposition 
of a random variable F G D 1,2 : Let H n the n-th Hermite polynomial, n G 
N, and denote by C n the closed linear subspace of L 2 (Q, J 7 , P) generated by 
the random variables {H n (B(Lp)),Lp G 7i, \\(p\\n = !}• Co will be the set of 
constants. Furthermore, denote by Q the a-algebra generated by the random 
variables B((p),ip G H. Then the space L 2 (£l,Q,P) can be decomposed into 
the infinite orthogonal sum of the subspaces C n : 

L 2 (n,g,p) = ®™ =0 c n . 

Moreover, denote by J\ : L 2 (Q) — > C\ the projection to the first chaos. If 
F G IB) 1 ' 2 such that sup s6 [ 0jl ] E| D S F\ 2 < oo, then we have 

J 1 (F) = 5(EDF) (12) 

almost surely, which is a straightforward consequence of the transfer principle 
for fractional Brownian motion, see e.g. Chapter 5.2 in [30] and Stroock's 
formula, see, e.g., Chapter 6 in [18]. 

2.2 Stochastic differential equations driven by fractional Brownian motion 

Throughout this article, we will impose the following assumptions: 
(Al) H> 1/2 
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(A2) a G C 2 (R), a G C 3 (R) with bounded derivatives, 
(A3) a, a bounded, inf^R a(x) > 0. 

The assumption (A3) is - a priori - required only for technical reasons. See 
Section 5.1 for a discussion. 

Under the assumptions (Al) and (A2) it is well known that 

X t = x + f a(X T ) dr + /* a(X T ) dB T , t G [0, 1], 

JO JO 

which is the integral equation corresponding to equation (1), has a unique 
solution X = (X t )te[o,i\ with 

E sup \X t \ p < oo 
te[o,i] 

for all p > 1. See e.g. [14], [27]. Here, the integral with respect to fractional 
Brownian motion is defined as a Riemann-Stieltjes integral. 

Moreover, we have X t G D 1,2 for all t G [0, 1] and 

D u X t = a(X u ) exp a'(X T ) dr + £ a'(X T ) dB^j l m (u), u, t G [0, 1]. 

(13) 

See [28], [31]. 

Since the diffusion coefficient is strictly positive due to the assumption (A3) 
we can use the Lamperti transformation, which will be an important tool 
throughout this article. Define 




Then d : R — > R is well defined, since a is strictly positive. Moreover, we have 




x G R. 



Note that d : R — > R is strictly monotone, thus the inverse function d : R — > R 
exists and satisfies 

(■&- 1 )'(x)=a(0- 1 (x)), xeR. 

A straightforward application of the change of variable formula for Riemann- 
Stieltjes integrals, see e.g. [37], yields that 

Y t = &(X t ), tG[0,l], (14) 
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is the unique solution of the stochastic differential equation 

dY t = g(Y t )dt + dB t , te[0,l], (15) 
Y = #(x ), 

with 

* w = ^Ur x e R ' 

Clearly, we also have 

X t = -&~\Y t ), te[0,l]. (16) 



Note that the mapping g : R — > R is twice continuously differentiable with 
bounded derivatives. 

Using (16) we can also give a different representation of the Malliavin deriva- 
tive of X t , t G [0, 1], which will be more appropriate for our purposes. See [12] 
for a similar representation in the case of Ito stochastic differential equations 
driven by Brownian motion. 

Proposition 1 We have 

D s X t = a(X t ) exp ^ (a' - (X T ) dr^j l m (s), s,t G [0, 1]. (17) 



Proof. Using the chain rule (9) we have 

D s X t = {d- l )'{Y t )D s Y t = a(X t )D s Y t , s,t G [0, 1]. 
Since (Y^) t6 [ 0i i] is the solution of equation (15), applying (13) yields that 

D s Y t = exp (j\'(Y T )d?j l m (s), s,t G [0,1]. 



Moreover, it holds 



g\x) = a'^\x))-—{d-\x)). 



Thus, the assertion now follows by (16). 

□ 
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In particular, the representation (17) implies that (D s Xt) Sj te[o,i] i s a bounded 
stochastic field. 



3 The degenerated case 



In this section we study under which conditions the solution of equation (1) 
is "degenerated" in the following sense: The solution X t at time t G [0, 1] 
does not depend on the whole sample path (B s ) se [o t t] of the driving fractional 
Brownian motion up to time t, but only on B t . As for Stratonovich stochastic 
differential equations driven by Brownian motion, see [36], this property can 
be completely characterized in terms of the drift- and diffusion coefficient. 

Theorem 2 There exists a mapping f G C 2 ([0, 1] x R; R) such that 

(X t )te[o,i] = (f(t,B t )) tem 
almost surely, if and only if (a'a — acr')(x) = for all x G R. 

Proof, (i) Suppose that a'a — aa' = 0. Using the Lamperti transformation, 
we have X t = ti' 1 {Y t ), t G [0, 1], with 

dY t = g(Y t ) dt + dB t , te[0,l], 
Y = #(x ) 



and 



Since 



a(v {x)) 



aa' . „ , , « . a'a — aa' 



a a 
the assumption a'a — aa' = implies that 

g'(x) = 0, xeR. 

Consequently, we have 

Y t = ^p\t + B t , te[0,l]. 
a(x ) 

Thus, we obtain the representation 

X t = f(t,B t ), t G [0,1], 
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with 



Due to the assumptions on a, we also have that / G C 2 ([0, 1] x R; R). 
(ii) Now assume that 

X t = f(t,B t ), te[0,l] (18) 

with / G C 2 ([0, 1] x R; R). From (18) and the chain rule for Riemann-Stieltjes 
integrals, see e.g. [37], we obtain 

X t = x + f f t (r, B T ) dr + f f x (r, B T ) dB T , t G [0, 1], 

JO JO 

almost surely. This yields 

f f t (r, B T ) - a(X T ) dr = - f f x (r, B T ) - a(X T ) dB T , t G [0, 1], (19) 
jo jo 

almost surely. Now set 

ZP = f ft(r, B T ) - a(X T ) dr, t G [0, 1], 

JO 

and 

Z? ] = f f x (r, B T ) - a(X T ) dB T , t G [0, 1]. 

JO 

Moreover, define the a-variation V"(Z) with stepsize 1/n of a stochastic pro- 
cess Z = (Z t ) te[0:1] by 



n-1 



V n( Z ) = J2 \ Z (i+l)/n - Z i/n 



i=0 



Clearly, we have by (19) that 

= C(^ (2) ) 

almost surely for all a > and n G N. Since the process (Z t ^) te [ 0) i] is pathwise 
continuously differentiable, we have that 

V l/H( Z (l^Prob Q 

as n — > oo. Moreover, Theorem 1 in [10] yields that 

V^(Z^) ^ E\B^ H f \Ur,B T )-a(X T )\^ H dT 

Jo 
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as n — > oo. Thus, we have that 

[ 1 \f x (r,B T )-a(X T )\ 1/H dr = 
Jo 

almost surely and it follows 

f x (t,Bt)=a(X t ), te[0,l], 
almost surely. This yields in turn that 

f t (t,B t ) = a(X t ), te[0,l], 

almost surely. Since X t = f(t, B t ), we can write the above two equations also 

as 

f x (t,B t )=a(f(t,B t )), te[0,l]. 

and 

f t (t,B t ) = a(f(t,B t )), te[0,l]. 
Since B t is a centered Gaussian random variable with variance t 2H , we obtain 

f x (t, x) = a(f(t, x)), f t (t, x) = a(f(t, x)), te (0, 1], x G R, 

which implies that 

O = f x t{t,x)-f tx (t,x)=ao , (f(t,x))-aa'(f(t,x)), t e (0,1], x G R. 

However, since a does not vanish due to our assumptions, the distribution 
of X t is absolutely continuous with respect to the Lebesgue measure with a 
strictly positive density for every t G (0, 1]. See e.g [28]. This implies that the 
image of the mapping / is the whole real line R. Now, we have finally that 

aa'(x) - aa'(x) =0, x eR. 

□ 

Thus, the solution of equation (1) is "degenerated" if and only if a and a 
commute in the usual sense of differential geometry. Since 

x[^ = E( Xi | Bi/ n , Bi/ n , . . . , B\ ), 

the following corollary is a straightforward consequence of Theorem 1. 

Corollary 3 Let (a' a — aa')(x) = for all xGl. Then we have 

(Epfx-X? ! 2 ) 1 / 2 ^ 

for all nGN. 
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Hence the above Corollary implies that X\ can be simulated perfectly - at 
least theoretically. The mapping /; [0, 1] x 1 -» 1 in Theorem 2 is given by 

f x (t, x) = a(f(t, x)), f t (t, x) = a(f(t, x)), t G (0, 1], 

with /(0, x) = Xo, x G R. The solution of this partial differential equation will 
be explicitly known only in some particular cases. 



4 The non-degenerate case. 



In this section, we determine the following lower and upper bound for the 

i n \ 

error of the conditional expectation X l in the non-commutative case: 
Theorem 4 It holds 

liminf ^^(ElXi-X^I 2 ) 172 > a H ^ \Ey t \ 2 dtj' 2 (20) 

and 

hmsup n H+l ' 2 (E|Xx -X^f) 1 / 2 < (3 H ( [' E\y t \ 2 dt)^ (21) 

where 

y t = (aa 1 - a'a)(X t ) exp Qf* a'(X s ) ds + j' a'(X s ) dB.^j , t G [0, 1] 
and «/f > and (3h > are too numerical constants depending only on H . 

Clearly, the random weight function y vanishes, if {a' a — aa') (x) = for all 
i6l. However, if equation (1) satisfies 



(nd) f 1 \Ey t \dt>o, 

Jo 



then the exact rate of convergence of the conditional expectation is n~ H ~ l l 2 , 
which is summarized in the following Corollary. Note that condition (ND) is 
satisfied, if e.g. (a'a — aa')(x ) ^ 0. 

Corollary 5 // (ND) holds, then there exist constants C\ = Ci(a, a, x , H) > 
and C 2 = C 2 (a, a,x ,H) > such that 

C x -n- H - x l 2 < (E|Xx -X[ n) | 2 ) 1/2 < C 2 -n-*- 1 ' 2 

for all n G N. 
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Consequently, the maximum rate of convergence, which can be obtained by 
an equidistant discretization of the driving fractional Brownian motion, is 
n~ H ~ l l 2 in this case. Moreover, we can now characterize the difficulty of equa- 
tion (1) in terms of its coefficients: If the drift- and diffusion coefficient com- 
mute, then Xi can be approximated perfectly, see Corollary 3. Otherwise, if 
a and a do not commute, then there are initial values xq G K. such that the 
exact convergence rate of the conditional expectation is n~ H ~ l l 2 . 

Theorem 4 fits in the known results for the case H = 1/2. In the case of a one- 
dimensional Stratonovich SDE driven by a Brownian motion W = (W t ) t6 [ 0i i], 
i.e. 

dV t = a(V t ) dt + a(V t ) dW t , V = v eR, te[0,l], 
it is well known that 

Urn „ • (E|V, - vi'W = j- (//E W «f .a)" 2 , (22) 
where in this case 

V} n) = E(V 1 \W i/n ,i = 0,...,n) 

and 

yf = (aa' - a'a)(V t ) exp Qf* a'(V T ) dr + £ a'(V T ) dW^j , t G [0, 1]. 

See, e.g. [7], [24] and [5]. In these articles, the key for the proof of (22) is to 
derive the asymptotic error distribution of n(Vi — v}^), which is a conditional 
normal distribution with zero mean and variance given in terms of a stochastic 
differential equation. Since this method essentially relies on the properties of 
Brownian motion, we could not imitate it and have to rely on indirect methods. 

4-1 Remarks 

(i) If the diffusion coefficient is constant, the boundedness of the drift coeffi- 
cient is not required in Theorem 2 and 4, since we do not have to apply the 
Lamperti transformation in this case. See also [22]. In particular, we obtain 
for the Langevin equation 

dX t = \X t dt + dB t , t G [0, 1], X = x 

that 

liminf n" +1 /2( E | Xl _x; n) | 2 ) 1/2 > a H W dtj 
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and 

1 /2 

limsup n H+1 ' 2 (E|Xx -X^ ] \ 2 fl 2 <Ph(C W dt) 

with 

y t = Aexp(A(l-t)), te[0,l]. 

(ii) Condition (ND) is a non- degeneracy condition for the first term of the 
chaos expansion of X x . Since (D s X 1 ) s6 [ 0)1 ] is bounded, we have by (12) that 

Xi = EX 1 + ( l ED s Xx dB s + ..., 
Jo 

If (ND) is not satisfied, then we have 

f 1 ED S X X dB s = c-B 1 
Jo 

for a constant c6i, which means that the first chaos of Xi is just a multiple 
of B 1 . Compare with Theorem 2, which yields, if a and a commute, the chaos 
expansion 

oo 

X ± = c + ci • B 1 + J2 d • HiiBi), 

i=2 

where q 6 R and Hi is the i-th Hermite polynomial. 

We strongly suppose that the conditions a'a—aa' ^ and (ND) are equivalent, 
i.e. either the error of the conditional expectation is zero or otherwise its exact 
rate of convergence is rT H ~ l l 2 . 

(iii) At first view, it may seem restrictive to consider only equidistant dis- 
cretizations of the driving fractional Brownian motion. However, for H ^ 1/2 
the increments of fractional Brownian motion are correlated and therefore the 
exact simulation of B tl , . . . , B tn is in general computationally very expensive. 
Given n iid standard normal random numbers, the Cholesky decomposition 
method, which is to our best knowledge the only known exact method for 
the non-equidistant simulation of fractional Brownian motion, requires still 
0(n 2 ) arithmetic operations after precomputation of the factorization of the 
covariance matrix. However, if the discretization is equidistant, i.e., tj = i/n, 
i — 1, . . . , n, the computational cost can be lowered considerably, making use 
of the stationarity of the increments of fractional Brownian motion. For ex- 
ample, the Davies-Harte algorithm for the equidistant simulation of fractional 
Brownian motion has computational cost 0(nlog(n)), see e.g. [9]. 

For a comprehensive survey of simulation methods for fractional Brownian 
motion we refer to [8]. 
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(iv) For multi-dimensional Stratonovich SDE, i.e., 

m 

dV t = a(V t )dt + J2v {i) (Vt)dW® 
i=i 

with a, : W 1 — > M. d , i = 1, . . . ,m and m independent Brownian motions 
i = 1, . . . ,m, the optimal rate of convergence, which can be obtained 
by point evaluation of the driving Brownian motions, depends on whether the 
diffusion coefficients commute or not. If they commute, that is, if we have 

E^^^E^^S U = l,...,m,/ = l,...,d, (23) 

k=1 UXk k=1 CXfc 

then the optimal rate of convergence is n~ l as in the one-dimensional case, 
where n is the number of evaluations of the driving multi-dimensional Brow- 
nian motion. (Here denotes the k-th component of o~^.) However, if (23) 
is not satisfied, then the optimal rate of convergence is rT 1 ! 2 . See e.g. [7]. 

The prototype example for the latter case is the two-dimensional Stratonovich 
SDE 

dV^ = V^dW^\ =0, 
dV t {2) = dW t {2) , \/ (2) =o. 

Clearly, we have = Jq dw} 1 ^ and it is well known that 

(E| V® - E( | W% W% i = 1, . . . , n )| 2 ) 1/2 = • n~^, 

see e.g. [7] 

In a forthcoming paper we will study, whether this phenomenon also appears 
for SDEs driven by fractional Brownian motion. 



4-2 Proof of the lower error bound in Theorem 4 



The Wiener chaos decomposition described in Subsection 2.1 will be the key 
for the proof of the lower error bound, since it provides a linearization of the 
problem, which we have to analyse. 



Proposition 6 We have 



E|Xi - X\ '\ z > E 



W|2 ^ ™ £ (B t - E(B t | B i/n ,i = 0, . . . ,n)) Ey t 



dt 



(24) 
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Proof. Clearly, we have that 

J 1 (X 1 )=5(EDX 1 ), 
see equation (12) in Subsection 2.1. Recall that 



D S X, = a(X^ exp ^ (a' - (X T ) di 



by Proposition 1 and define the mapping m : [0, 1] — > ffi. by 

m(s) — ED S X ± , se [0,1]. 



Since a and a' — ^- are bounded due to our assumptions, it follows by dom- 
inated convergence that the mapping m : [0, 1] — > R is continuously differen- 
tiate with 

m'(s) = EDA| fl '~|(X ( ), se[0,l]. 
Using (13) we obtain 

m'(s) = Ey s , se[0,l]. 

Since the Skorohod integral and the Riemann-Stieltjes integral coincide for 
smooth deterministic integrands, see (10), we finally have that 



J 1 (X 1 ) = f 1 E D t X 1 dB t = B 1 Ea(X 1 ) - C B t Ey t 
Jo Jo 



dt. 



The projection of X l =E(X\ \ Bi/ n , i — 0, . . . , n ) to the first chaos is given 
by the second term of the Hermite series expansion of the conditional expec- 
tation: 

n 

X 1 = EX ± + ^ a iBi/n + • • • 
i=l 

with Oj e M. for i — 1, . . . , n. See e.g. [18]. Hence we have 

J l (X i ? ) ) = Y j a l B l/n 
i=i 

and it follows 



E|Xi -xS n) | 2 > E 



BiE<t(Xi) - 5 t E^ dt-J2 a l B l/n 
Jo i=i 
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by linearity of the projection to the first chaos. The term on the right hand 
side is the error of the quadrature formula 

n 

T n (B) = B 1 Ea(X 1 )-Y / a l B l/n 

i=i 

for the approximation of the integral 

i(B) = f 1 B t Ey t dt. 
Jo 

Clearly, by definition of the conditional expectation the best quadrature for- 
mula, which uses B , B 1 / n , . . . , Bi, is given by 

l n {B) = f 1 E(B t | B i/n , i = 0,...,n)Ey t dt. 
Jo 

Thus, we finally obtain 



E\Xx -X[ n) | 2 > E 



f 1 B t Ey t dt - f 1 E(B t | B l/n , i = 0, . . . , n) Ey t dt 
Jo Jo 



□ 



To finish the proof of the lower error bound in Theorem 4, we now have to 
analyse the quantity 

e(n) = (e £ Ey s (B s - E(B s \B l/n , i = 0, . . . , n))ds j . 

For this, we recall some well known facts about reproducing Kernel Hilbert 
spaces and weighted integration problems for Gaussian processes. See e.g. [32] 
and the references therein. 

Let W = (Wt)tgK be a continuous stochastic process with covariance function 
(K(s, i)) s ,teR. Then the reproducing Kernel Hilbert space (Sj(K), (•, -)k) cor- 
responding to the process W is the uniquely determined Hilbert space of real 
valued functions on R such that 

K(.,t)e?)(K) 

and 

(h,K(-,t)) K = h(t), teR 

holds for all h e f>(K). 

Now consider the linear functional 

I(W) = I p(t)W t dt, 
Jr 
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where the function p G C(R) is non-negative and has compact support. Then 
the error of a quadrature formula 

n 

T(W) = ]T a, W tl 

i=i 

with aj, tj G H. for « = 1, . . . , n, can be characterized in terms of the reproduc- 
ing Kernel Hilbert space as follows: 



(E|/(W)-/(W)| 2 ) 1/2 = sup 



p(t)h(t)dt-J2^h(t t ) 
i=i 



;heSj(K): \\h\\ K <l 



In addition, if the process W is Gaussian and if we use the conditional expec- 
tation of the functional given the evaluations of the random process, i.e. 

7(W) = E(/(W)|W^ = l,...,n), 

as quadrature formula, then we have 

(E|/(W) -7(W) 2 ) 1/2 

= sup | J R P(t)h(t) dt ; h G S)(K) : \\h\\ K < 1, h(U) = 0, i = 1, . . . , nj . 

See e.g. section III. 2 in [32]. Now we can finish the proof of (20): 

Proof of the lower bound in Theorem 4: By the above considerations it 
remains to analyse the quantity 

e(n) = sup| J\(t)h(t)dt ; h e Sj(R H ) ■ \\h\\ RH < 1, h(i/n) = 0, i = 0, . . . , nj 

with p(t) = E^i, t G [0,1]. Here ($)(Rh), ('■>') Rh) ^ s ^ e reproducing kernel 
Hilbert space of B = (B t ) teR , which is the fractional Brownian motion defined 
on the whole real line. 

It is well known that 

G = {g G C°°(R) : suppg compact, ^ suppg} 

is dense in $)(Rh) and that 

\\h\\l H = c H f \ur^\Hh)(.u)\ 2 du 

for h G G, where the constant ch > is known explicitly, see e.g. section 6.1 
in [32] and [33]. Here F{h) denotes the Fourier transform of h, defined by 



!F(h)(t) — / exp(— itu)h(u) du. 



18 



Now consider the stationary Gaussian process V = (Vt)tm with covariance 
kernel 

K(0,t) = [ exp(itu)f(u)du, tel, 
where the spectral density / : H. — > R is given by 



/( " ) = 4^ (1 + " 2) """ I/2 ' " £ 



(25) 



Then the reproducing Kernel Hilbert space corresponding to V is the Bessel 
potential space 



f8 H+1 ' 2 = {heL 2 (R):\h\ H+1/2 <oo} 



with 



\h\% 



4tt 2 



f\uy l \F(u)\ 2 du. 



{he<S H : h(0) = 0} C $)(R H ), 



We have 
and 

\\h\\ RH < \h\ H 

for h £ 23 see e.g. section 6.1 in [32]. Hence it follows 



e(n) > sup | j\(t)h(t)dt ;he*B H : \\h\\ R[I < 1, /i(i/n) = 0, i = 0, . . . , nj 



and moreover 



e(n) > sup | J\{t)h{t)dt ;/ie*8 F : < 1, /i(i/ra) = 0, i = 0, . . . , nj . 



However, the quantity on the right hand side of the above equation corresponds 

t0 1/2 

e(n)= ^e\J* p(t)(V t -E(V t \V i/n ,i = 0,...,n))dt\ ) , 

where V is the stationary Gaussian process with covariance kernel given by 
(25). 

From [34], Proposition 2.1, we now obtain that 



hm «pl = 1 

™- +0 ° s(n) 



with 



s(n) 



£ f(u + n2nj)\F(p(u))\ 2 du. 

-717T rr,l 



jez\{o} 
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We have 

/nir 
\F{p{u))\ 2 du- ]T f^ + n2nj) 
jez\{ } 

with £ G (— nn,mr). Note that by symmetry of the spectral density it holds 

oo 

f(u + n2nj)\ = 2j2f(u + n2nj), ueR. 
jez\{o} j=i 

Since / is strictly decreasing, it follows 

00 run 

s(n)>2^/K27ri + 7r)) / \F{p{u))\ 2 du. 

7=1 J-nir 



Now consider the quantity n 2H+1 s(n). Clearly, we have 

liminf n 2/m s(n) > — — liminf V (n- 2 + n(2j + if)-*- 1 / 2 / \F(p(u))\ 2 du. 
n^oo v 1 - 2ir 2 c H ^ v 7m 1 v v 

Fatou's Lemma and Parseval's equality yield 

liminf n 2H+1 s(n) > — V- — — — r^ry— 7 / |pH| 2 *x. 
n-^> v y - c H 2n 2 ^ ((2j + l)^) 2 ^ 1 Jr lHy n 



Since 

liminf n 2H+1 e(n) > liminf n 2H+l s(n), 

it follows 



liminf n 2/m e(n) > c& / \p(u)\ 2 du 

n— >oc Jjg 



with 

1 1 ^ 1 

^ 27i 2H + 2 c H f^ (2j + l) 2 ^ 1 ' 

Since p(fj = E^, t G [0, 1], and 

e(n) <E|Xx-xS n) | 2 
by Proposition 6, we finally have shown the lower bound in Theorem 4. 

□ 

4-3 Proof of the upper bound in Theorem 4 

For the proof of (21) we will again use an indirect method. 
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Denote by = (B^ n) ) tem the piecewise linear interpolation of B based on 
an equidistant discretization with stepsize 1/n, that is 

B ( t n) = B k/n + (nt - k) (S( fc+ i)/ n - fife/™) , t G [k/n, (k + l)/n], 
and consider the following stochastic differential equation 

dx\ n) = a(xi n) ) dt + a(xi n) ) dB[ n \ t G [0, 1], 

V"( n ) _ „, 

A — Xq. 

This equation has a pathwise unique solution, since B^ has finite total vari- 
ation for every n G N. In particular, we can again apply the Lamperti trans- 
formation given in Subsection 2.2 and obtain that 

5J») = ir i(yW) > tG [ ,i], 

where Y^> = (y^) te [ 0;1 ] is the unique solution of 

dY t {n) =g(Y t {n) ) dt + dBl n \ te[Q,l], 
Yt ] =t?(x ). 



Clearly, we have 



El^-XS'V^ElXi-X^I 2 



for every n G N. The next Proposition characterizes the leading term of the 
error of the approximation x[ : 



Proposition 7 We have 



exi - xr\ = E 



2 

0)|2 



f 1 y t (B t - Bl n) ) dt 

Jo 



+ 0(n- 21i ). (26) 



Proof. We will denote constants, which depend only on xo, H, a, a and their 
derivatives by c, regardless of their value. 

(i) We first we establish the following estimate: 



sup \Y t \+ sup \Y t (n) \ < c sup \B t \. (27) 
te[o,i] te[o,i] *e[o,i] 



For this note that 

r* 
'o 



^| < \^ 1 (x )\ + f g(Y T )dr+\B t \, t G [0,1]. 
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Since g is continuously differentiable with bounded derivative, it satisfies a 
linear growth condition and we obtain 

sup \Y S \ <c + c ( sup |y s |dr+ sup \B S \, te [0,1]. 
se[o,t] Jo se[o,r] se[o,t] 

Thus, Gronwall's Lemma yields that 

sup \Y t \ < c sup \B t \. 

*6[0,1] *6[0,1] 

Similar, we obtain 

sup \Y t {n) \ < c sup \B t \, 

te[o,i] te[o,i] 

since 

sup |.B^| < sup \B t \. 

te[o,i] te[o,i] 

(ii) Using the Lamperti transformation, we can write 

x t - x, (ri) = $~\Y t ) - r^rt, f G [0, 1]. 

Since we have 

(iT 1 )'^) = a(^- 1 (a;)), 
= (7( 7 / (^- 1 (a;)), 

we obtain 

X t - X[ n) = a(X t )(Y t - Y t {n) ) + l -aa'^-\e t + (1 - 6 t )Y< n) ))(Y t - f t (n) ) 2 

for t G [0,1] with a random 9 t G (0,1). So we only need to consider the 
difference between Y and Y^ n \ For 

Z t = Y t -Y t (n \ t G [0,1], 

we obtain 

z t = Bt - £j n) + /' ^(y T ) - j(y T W) dr, f G [0, 1]. 
The Lipschitz continuity of (/ implies that 

\Z t \ < \B t - B[ n) \ + c /* |Z T | dr, t G [0, 1] 

Jo 

and Gronwall's Lemma yields 

\Z t \ < c f \B T -B^\dT, te[0,ll. 
Jo 
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Thus, it follows 



sup E\Z t \ p <c-n- Hp , (28) 
te[o,i] 



since clearly 

sup E\B t -Bi n) \ p <c-n- Hp . 

*e[o,i] 

Using the boundedness of a and a' we obtain 

X t - Xl n) = a(X t )Z t + itf (29) 

with 

te[o,i] 



sup E|i2t (1) | p < c-n" 2Hp . 



(hi) Now, we analyse the process Z = (Z t ) t6 [ 0i i] in more detail. We can write 

Z t = B t - B\ n) + /' g'(Y T )Z T dr + R?\ t G [0, 1], (30) 
Jo 

where 

R? ] = [\g'(0 T Y T + (l - fl T )y r W) - ^(y T ))z T dr, * g [o, 1], 

with # r G (0, 1). Since is Lipschitz continuous, it follows 

g'(6 T Y T + (1 - 9) T Y^) - g\Y T )\ < c\Z T \, r G [0, 1], 
and consequently 

\R? ] \ < c / \Z T \ 2 dr, t G [0,1], 
Jo 

which in turn yields together with (28) that 

sup E|i4 2) r < c-n~ 2Hp . (31) 

te[o,i] 

Now consider the equation 

Z t = B t - B[ n) + /* g\Y T )Z T dr, t G [0, 1]. 
Jo 

Its unique solution is given by 

Z t = j\'( Y r) exp (/V(^) ds) (S T - flW) rfr, t G [0, 1]. (32) 
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Since 



Z t -Z t = R? ] + f g'(Y T )(Z T - Z T ) dr, 
Jo 

again an application of Gronwall's Lemma and (31) yield that 



sup E\Z t - Z t \ p < c ■ rT 2Hv . (33) 
te[o,i] 

Therefore, we obtain from (32) and (33) that 

Z t = f Q 9'(Y T ) exp ( jf g'(Y s ) ds^j (B T - B™) dr + itf 

with 

sup E\Ei 3) \ p <c-n~ 2Hp . 
te[o,i] 

(iv) Using (29) we finally have that 

X t -X[ n) = a(X t )Z t + R?\ fe[0,l], 

with 



sup E\Ri A) \ p <c-n- 2Hp . 

te[o,i] 



Since 



aa' 



g'{Y s )= («'- — ) se [0,1], 



we obtain 

ft / nrr' 



(4) 



x exp ^jf* ^a' - ^ ds^J (S T - 5| n) ) rfr + i? 

for t £ [0,1], and the assertion follows, since 

a(X t ) exp ( jf* (a' - ^) (X s ) <fe) l [0 , t] ( S ) 

= a(X r ) exp QT a'(X s ) rfs + J' a'{X s ) dB^j l m {s) 

by (17). 

□ 
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Now, it remains to analyse the quantity 



E 



ty t {B t -Bl n) ) 

Jo 



dt 



2 



Note that this is the mean square error of the trapezoidal-type quadrature 
formula 

1(B) =n Pfa-tWtdt-Bo 
Jo 

71—1 / rt i+ i rti \ 

+ n E / 1 - Wt dt+ 1 (t- U-jyt dt ■ B l/n 

i— i \Ju Ju-i J 

+ n f 1 (t-tn-Mdt-Bu 
Jt„-i 

where U — i/n, i — 0, . . . , n, for the approximation of the weighted integral 

i(B) = Cy t B t dt 

Jo 

with random weight function y = (3^)te[o,i]- 

Combining the next Proposition and Proposition 7 we obtain the upper bound 
of Theorem 4. 

Proposition 8 We have 



lim n 2H+1 E 



f 1 y t (B t - B[ n) ) ds 2 = \C(-2H)\ f 1 E\y t \ 2 dt, 

Jo Jo 



where ( denotes the Riemann Zeta function. 

The proof of this Proposition is rather technical. Therefore we postpone it to 
the Appendix. Note that for the constant (3h in Theorem 4 we have 

/% = |C(-2#)|. 



5 McShane's scheme method 



Theorem 4 states in particular that the rate of convergence of the conditional 
expectation is at least n~ H ~ l l 2 . Clearly, the conditional expectation is ex- 
plicitly known only in some exceptional cases and thus is not an implementable 
approximation scheme in general . In this section, we present a feasible ap- 
proximation scheme, which is almost as good as the conditional expectation 
in the sense that its convergence rate is also at least n~ H ~ l l 2 . 
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The McShane's scheme for stochastic differential equation driven by fractional 
Brownian motion is defined by = x and 

Xl% = + a(*J°)A + a(X^)A k B + \aa\xP){A k B) 2 (34) 

+ \{oa> + a'a)(X^)A k BA + l -aa'(X^)A 2 

+ ^(aV + ff (,') 2 )(f ) )(A fc B) 3 

for k = 0, ...,n — 1, where A = 1/n and A k B = B( k+1 y n — B k / n . For 
Stratonovich SDEs driven by Brownian motion, this scheme was studied e.g. in 
[17], [24] and shown to be asymptotically efficient, i.e. 

lim nE|fA -Z7S n) | 2 = lim nE^ - U^\ 2 = — C ny^f dt. 

n-*oo 1 n-»oo 12 Jo 

The following Theorem gives in the non-degenerated case the exact conver- 
gence rate for McShane's method for SDEs driven by fractional Brownian 
motion: 

Theorem 9 For the approximation scheme given by (34) we have 

Jim n^/^ElXx- X^\ 2 ) 1 ' 2 = Ph (j dt) . (35) 

Note that the asymptotic constant on the right hand side of (35) vanishes, if 
and only if (a! a — aa')(x) = for all x G K. See Remark (iii) in Section 5.1. 
Hence, if condition (ND) holds, then both McShane's method and the condi- 
tional expectation have exact rate of convergence n~ H ~ l l 2 and thus McShane's 
method is optimal in this case. 

In the degenerated case, we obtain the following upper bound for the error of 
McShane's scheme. 

Proposition 10 // (a' a — aa')(x) = for all x E R, then there exists a 
constant C = C(a, a, x , H) > such that 

(E|Xx- X^\ 2 ) 1 ' 2 <C-n~ 2H 

for all nGN. 

Proof. This is a straightforward consequence of Proposition 7. 

□ 



26 



Thus, in this case McShane's scheme has a convergence rate of at least n 2H . 
For the case of a zero drift coefficient, i.e. a = 0, it is shown in [20] that 



n 



4H-1 



Prob 



—aiX,) j (V :5 \aa'a" - aV")( X,) <ls 



for n — > oo. So we strongly suppose that the exact convergence rate of Mc- 
Shane's scheme in the degenerated case is n~ 4H+1 . 



5. 1 Remarks 



(i) In this article, we use indirect methods to determine the mean-square 
error of the considered approximation schemes. The main reason for this is 
that moment estimates for Riemann-Stieltjes integrals driven by fractional 
Brownian motion are much more involved than for Ito integrals with respect 
to Brownian motion. 

Recall that, if the process u = (u t )te[o,i] satisfies appropriate smoothness con- 
ditions, the relation between the Riemann-Stieltjes integral and the Skorohod 
integral is given by 

f 1 u t dB t = 5{u) +a H f 1 C D s u t \s - t\ 2H ~ 2 dsdt (36) 

JO JO JO 

and moreover 

E\5{u)\ 2 = a H I Eu s u r \s -r\ 2H - 2 drds (37) 
^[o,i] 2 

+ a\ ( ED r u s D r ,uJr - s'\ 2H ~ 2 \r' - s\ 2H - 2 drdr'dsds' . 
Ao,i] 4 



Since in both expressions the Malliavin derivative appears, it is not possible 
to use them for a direct error analysis. To illustrate this, consider e.g. the 
continuous Euler method for equation (1), which is given by X§ = x and 

X? = X* n + f a(X* n ) dt + f a(X* n ) dB t , t G [k/n, (k + 1)/N) 

Jk/n Jk/n 

for k — 0, . . . , n — 1. Here we have 

X t - Xf = jf* a(X t ) - a(Xf nt]/n ) dt + J* a(X t ) - a(Xf nt]/n ) dB t , t G [0, 1]. 

Applying (36) and (37) would yield an equation, which involves the first Malli- 
avin derivative of X t — Xf for t G [0, 1]. Thus, for analyzing the Euler method 
in this way, we would need to control the difference between the Malliavin 
derivative of the solution and the Malliavin derivative of the Euler method. 
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But this involves the second Malliavin derivative etc. and we cannot have 
closable formulas. 



In this article, we apply the Lamperti transformation to avoid this problem. 
Essentially the Lamperti transformation reduces the error analysis of the con- 
sidered equation to the error analysis of a related equation with constant 
diffusion coefficient, for which the above problem does not appear. The price 
we have to pay for this procedure is the quite strong assumption (A3) on the 
drift- and diffusion coefficient, which ensures the integrability of the remainder 
terms in the error analysis. 

If the pathwise error is considered instead of the mean square error, then 
assumption (A3) can be avoided. Under the assumptions (Al) and (A2) and 
using the Doss-Sussmann transformation as reduction method to analyse the 
Euler scheme, it is shown in [23] that 



n 2H ~ 



- X? ) ^ -\ C o'(X B )D s X 1 ds 
2 jo 

for n — > oo. 

(ii) For the proof of Theorem 9 we will use Proposition 7, i.e. 



E\X 1 -Xl n) \ 2 = E 



f 1 y t (B t - B[ n) ) dt 

Jo 



2 + 0(n~ 2H ), 



where 



dx\ n) = a(x| n) ) dt + a(xi n) ) dB[ n) , t € [0, 1], (38) 

and is the piecewise linear interpolation of B with stepsize 1/n. Therefore, 
it remains to compare x[ n ^ and X^ to show Theorem 9. Note that equation 
(38) actually corresponds to a system of piecewise random ordinary differential 
equations: If we define 

x {k \t) = a(x^(t)) + n(B {k+1)/n - B k/n )a(x^(t)), t e (k/n, (k + l)/n] 
x^ k \k/n) =x i - k ~ l \k/n) 

for k — 0, ... n — 1 with x^(0) = xo, then we have 

Xl n) =x i - k \t), t G [k/n, (k + l)/n]. 

Hence Theorem 9 will also hold for any method, which approximates the 
above system of piecewise ordinary differential equations with convergence 
rate n" H ~ 1 / 2 ~ e for e > 0. Compare with [6] for the case H = 1/2. 
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(iii) Since 

y t = {aa> - a'a){X t ) exp Qf* a'{X s ) ds + jf* a'(X s ) dB^j , t G [0, 1], 

the condition (acr' — a'<r )(x) = for all x G R clearly implies j*^ E|^ t | 2 = 0. 

On the other hand, assume that Jq 1 E|X| 2 dt = 0. However, Proposition 1 now 
yields that 

/ E\(aa' -a'a)(X t )\ 2 dt = 0. 
Jo 

Since a is strictly positive, the distribution of X t is absolutely continuous 
with respect to the Lebesgue measure with a strictly positive density for every 
t G (0, 1]. See e.g [28]. Hence it follows (aa' - a'a)(x) = for all xeR. 

Consequently, we have Jq E|3^| 2 dt = if and only if (aa' — a'a)(x) = for all 
x G R. 



5.2 Proof of the convergence rate of McShane's method in the non- degenerated 
case 



To show (35), we require the following Lemmata, which both can be shown 
by straightforward calculations. The first one gives a Taylor expansion of the 
Lamperti transformation, while the second considers the regularity of the drift 
coefficient of the reduced equation: 

Lemma 11 For the mapping $ : R — > R given by 

r x i 



#(z) = r -ft\dt, 

Jo ait) 



we /iawe i? G C 4 (R; R) wii/i bounded derivatives. In particular, it holds 



a(x) ' 



$"(x) = 



a'{x) 
a 2 (x) ' 



o" 2 (a;) a" 3 (re) 



for x G 



Lemma 12 Define g : R — > R fry 



<7 
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Then g is bounded and we have g G C 2 (R; M.) with bounded derivatives. More- 
over, it holds 

Proof of Theorem 9. We will denote constants, which depend only on xq, 
H, a, a and their derivatives by c, regardless of their value. Using Proposition 
7, the assertion of Theorem 9 follows, if we show that 

E\x[ n) -X[ n) \ 2 <c-A 4H . (39) 

(i) For this, define 

Y k = #(xl n) ), k = 0,...,n, 

which turns out to be an approximation scheme for Y^: A Taylor expansion 
yields 

+ V(x^)(xf + \-x^) 3 

+ + (1 _ e k )xit)(xl% - x ( k n) r 

for k — 0, . . . , n — 1 with a random 6 k e (0, 1) . Since a and a are bounded 
together with their derivatives, we moreover have that 

\XI% - Xi n) \ < c • (A + \A k B\ + \A k B\ 2 + | A k Bf) . (40) 

Hence we obtain that 

n+i = Y k + #\X^)(xi% - Xl n) ) (41) 

+ V(x^)(xa-xfv 
+ V(x^)(xa-x^) 3 +4 i) 

for k — 0, . . . , n — 1 with 

sup E\R [ k 1] \ p < c- A 4pH 

k=0,...,n— 1 

by the boundedness of and the estimate (40). 
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Now we have to analyse the different terms of the above recursion scheme. For 
the first term we have 

tf'(^ n) )(*a - = ^ n) )A + A k B + l -a'{X^){A k Bf (42) 



+ 1 -(**" + (*>) 2 )(xi n) )(A k B)\ 



while for the second term we obtain 

lnxl n) )(xl% - = -^(^)A 2 - (A,5) 2 (43) 

- ^(xl n) )A k BA - \{a'f(X^){A k Bf + i#> 

<J A 

with 

sup E\R k 2) \ p <c-A^ 1+2H \ 

fe=0,...,n-l 

Finally, for the third term we have 

l^'(xi n) )(n% - XP? = ~\ - 2(^) 2 ) (X^)(A k B)t + R k 3 \ (44) 
where 

sup E\R k 3) \ p <c-A^ 1+2H \ 

k=0,...,n— 1 

Combining (42), (43) and (44), it follows 
with 



sup E|4 4) |p <c-A^ 1+2H ). 

fc=0,...,n— 1 

Using Lemma 12 we can write the above recursion as 

Y k+1 = Y k + A k B + g{Y k )A + ^'(y fc )A fc SA + \gg\Y k )A 2 + 4 4) (45) 
with k — 0, . . . , n — 1. 

(ii) Now it remains to analyse the difference between y( n ) and F. For this set 
U — iA for i = 0, . . . , n. Moreover, recall that is given by 

? ) ^(? ) ), f e[0,i], 
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and that satisfies the integral equation 



Y t (n) = + jf g(y}») ds + A t B^ t g [U, ti+i], 



with Yq H ^ = $(xo). A Taylor expansion yields 



(46) 



for fc = 0, . . . , n — 1 with 



/ j T g{g'g)'{Y^)dudrdt 

^k ^k 

Jt k Jt k Jt k A 

+ p ,« rM } . j|f? ) )Writ 

Jt fe Jt fe Jt fe A z 

Straightforward computations and the boundedness of g, g' and g" yield 

sup E|4 5) |f <c-A^ 1+2ff ). 

fc=0,...,n— 1 



(3) Now set e k = Y^ — Y k for k — 0, . . . , n. Using (45) and (46) we have 



e k+1 = e k + (g(Y t ^)-g(Y k ))A 

l(g'(Y t M)-g'(Y k ))A k BA+ 1 - 



+ L -{g\Y t ^) - g'{Y k ))A k BA + -{gg'{Y^) - gg'(Y k ))A 2 + itf 



with e = and 



sup E\R k 6) \ p <c- A p ( 1+2H \ (47) 

k=0,...,n— 1 



Since g, and gg' are Lipschitz continuous due to Lemma 12, it follows that 

|e fc+ i| < \e k \(l + cA + c sup |B t |A|) + \R k 6) \ 

te[o,t] 

for /c = 0, . . . , n — 1. The discrete version of Gronwall's Lemma yields 

n— 1 

|e„| <exp(c(l+ sup |5 t |)) £ |4 6) I 
te[o,t] fe=0 



32 



and (47) implies that 

E|e n | p < c- A 2pH . 

Now (39) follows from 

\x[ n) - xi n) \ = |y^(n (n) ) - <P~\Y n )\ < c • \e n \. 



□ 



A Appendix 

In this section, we show Proposition 8 following section 3.5.3 in [22], where 
the case a = const, is studied. 

Recall that we have to show 
lim n 2H+1 E 



['y^Bt-B^ds ^\a-2H)\f\\y t 

Jo Jo 



2 dt. 



For this, we will need the Malliavin derivative of y = (D^)t e [o,i] and some 
technical Lemmata. Note that the weight function y can be written as 

y t = a{X x ) (a* - (X t ) exp ^ (a' - (X T ) dr^, te [0, 1], 

see Proposition 1, and that 

D s X t = a(X t ) exp ^ (a' - (X T ) dr^j l m (s), s, t e [0, 1]. 



Moreover, denote 



h(x) = la' j (x), x e 



for notational simplicity. 

The next Lemma follows by straightforward computations and the assump- 
tions on a and a. 



Lemma 13 For all p > 1 there exists a constant K{p) > such that 

E\y t -y s \ p <K( P )\t-s\p H 

for all s,t e [0, 1]. 
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The following Lemma can be shown by using the product rule (8) and the 
chain rule (9) in Section 2.1 for the Malliavin derivative. 

Lemma 14 We have y t G D 1 ' 2 for all t G [0, 1] with 

D s y t = (^j (x l )y t D s x l + h'(x t )D a x t D s x 1 + y t J' ti{X T )D s X T dr 

for s,t G [0, 1]. In particular, there exists a constant K > such that 

sup \D s y t \ < K. 
s,te[o,i] 

Denote in the following A = 1/n and 



U = iA, 



i — 0, . . . , n, 



U+i/2 = -(U + t i+1 ), i = 0, ...,n-\. 

We will also use the notation instead of y ti and instead of D s y t . in 

what follows. Moreover, set 



</>{s,t) = H(2H - 1) • \s - t 



2H-2 



s,t G [0,1]. 



Finally, we will again denote constants, which depend only on xq, H and a, a 
and their derivatives by c, regardless of their value. 



Lemma 15 We have 



sup 

i,j=0,...,n-l 



jf 3+1 J q D s yi ■ (t - t j+1/2 ) ■ 0(s, t) ds 



dt 



< c-A 2H+1 . 



Proof. We have 

rti+i rl 



J t 3+1 J o D s y t ■ (t - t j+1/2 ) ■ 0(s, t) ds dt 



D s yi ■ (t - t j+ i/ 2 ) ■ 0(s, t) ds dt 

tj+l ftj+2 



rt j+1 ftj+2 

+ / / D s y t ■ (t - t j+1/2 ) ■ 0(s, t) ds dt 

J tj j — 1 

+ f J+1 f D s y t ■ (t - t j+1/2 ) ■ <j)(s, t) ds dt, 

Jtj Jtj + 2 

with the convention that t_i = and t n+ i = 1. 



(A.l) 
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(i) We start by considering the second integral. Here we obtain 

rtj+i rtj+2 



f J+1 / 1+2 D s y, t ■ (t - t j+1/2 ) ■ <f)( s , t) ds dt 

Jtj Jtj-1 

rt j+1 rt j+ 2 

A- / <j>(s,t)ds 

Jtj Jtj — i 



< c 

using Lemma 14. Since moreover 

/ / </>(s,t)dsdt< c- A 2H , 

3 — 1 

by straightforward calculations, it follows 



dt 



f 3+1 f 3+2 D s yi ■ (t - t j+1/2 ) ■ <f)( s , t) ds dt 

Jtj Jtj — I 



<c-A 2H+1 . 



(ii) Now we study the first integral of the right hand side of (A.l). This integral 
clearly vanishes, if j — or j = 1. Now consider j — 2, . . . , n — 1. In this case 
< s < tj-i <tj<t< tj + i and thus we have 



f 1+ \t-t ]+1/2 )-(f)(t,s) dt<0 
Jtj 

for all s G [0,t J+ i]. The mean value theorem now implies that 
J t h+1 D s y t ■ (t - t j+1/2 ) • t) ds dt 



-3-1 



U 



j + 1 



D s yi (t-t j+1/2 )-(f>(t,s)dtds 
o Ju 



rtj+i rtj-i 

= ■ / (t - tj+1/2) ■ t) ds dt, 

J tj J 

where fj,ij is random and satisfies 



inf D s yi < fXij < sup D s yi 
s e[°.*i-i] se[o,tj_i] 



By partial integration we obtain 



-j+i 



(t-t j+1/2 ) -(f)(t, s)dt 



A ■ ((t j+1 - s) 2H ^ + ( tj - s) 2 ^) - H f t j+ \t - s) 2H ^ dt 



H 

f j t tj+ \t J+1 - sf H -^ + ( tj - *r-i - 2( t - s)™-* dt. 
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Integrating with respect to the variable s now yields 

(t — t j+ i/ 2 ) • (j)(s,t)dsdt 

= —— ({tj+i - tj_i) 2H - tj+i + (tj - £,-i) 2if - t 2H ^ 
r 1 ' 2H dt. 



Note that the above expression is the error of the trapezoidal approximation 
for the integral 



-3 + 1 



9j(t) dt 



with 



9j ■ [tj: t j+1 ] - R, 9j (t) = l'((t- t 3 ^f H - t 2H ) , t e [tj, t j+1 }. 



Since Qj G C 2 ([tj, tj+i\), it clearly holds 

rtj+i ftj-i 



■dt 



9'm _ A a 



12 



/■tj+i rtj-i 

J t J Q (t-t j+1/2 ) ■ <j>(s,t)dsi 

with £j G (tj,tj+i). Since 

g"{t) = 2H(2H - 1) ■ ((t - t^) 2H - 2 - t 2H ~ 2 ) , t G fe,f,- + i), 



we obtain 

So finally, it follows 

rtj+i ftj-i 



sup \g'j(ij)\ < c • A 



2.H-2 



j=2,...,n-l 



A^i • (t - ^+1/2) • (it 



<c-A 1+2H . 



(iii) It remains to consider the third integral, which vanishes for j = n — 2 and 
j — n — 1. Here we have < < t < tj +1 < tj +2 < s < 1 and thus 



f tj+ \t-t j+1/2 )-(f)(s,t) dt>0 



for all s G [i,-+2, 1]- Hence we get by the mean value theorem for integration 

D s yi ■ (t - £7+1/2) • 0(s, t) ds dt 
t j+ i ri 



tj J tj+2 



fiij ■ I I (t- t j+1/2 ) ■ 0(s, t) ds dt, 

Itj Jtj+2 
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with 



Now, the term 



inf D s yi < fiij < sup D s y(ti). 

se[tj+2,l] S£[t j+2 ,l] 



/ / (t-t j+1/2 )-<j>(s,t) 



Itj Jt j + 2 

can be treated analogous to (ii). 



dsdt 



□ 



Define 

Wj (t) = t - tj+i/2, t E [tj,t j+1 ], 

for j = 0, ... ,n — 1. In the next Lemma, we precompute the error of the 
weighted integration problem. 



Lemma 16 We have 







1 _ 2 \ 1/2 
B t -B^)y t dt 



n—1 n—1 



i=0 j=0 



= ( E E E y&i [ U+1 f t 3+1 Wj{t) Wi {s)<j>{s, t) dtds} + 0(A 2H ). 



1/2 



Proof. First note that 

u 



2\ 1/2 



f\ Bs -B^)y s ds 

Jo 

n—1 ft+i ~ 
„-_n 



2\ 1/2 



+ 0(A iW ) 



by Lemma 13. Moreover, applying partial integration we get 



rtj+i 

Jti 



B t - B[ n) dt 



■j+i 



Wj (t) dB t 



for j = 0, . . . , n — 1. By relation (10) in Subsection 2.1 and Lemma 14 it follows 

y j [ tj+1 Wj (t)dB t = f 1 D^jwjit^is^dsdt + Siyjwjl^.^). 



By Lemma 15 we have 

rtj+i fl 



1 3+1 I D s y j w j (t)(f)(s,t)dsdt 
Jti Jo 



<c-A 2H+1 . 
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Hence, it follows 



E 



=o 



2\ 1/2 



n-l 



E 



2\ 1/2 



+ 0(A 



2H\ 



Using the isometry (11) for Skorohod integrals with respect to fractional Brow- 
nian motion, we obtain 



e (s (y iWi l [ti , ti+l] ) 5 (y jWj l [t . >t 



E 



/ / 3^ • Wi^wjit) ■ 0(s, t) dt ds 
+ E/ / / / ^^(ti)-^^-^) 

Jti Jtj j0 JO 

x 0(si, t 2 )0(s 2 , ti) ds 2 rfsi dt 2 dti 



for i, j = 0, . . . , n — 1. Since fortunately 

rtj+i ftj+i rl r~t 



Ds^iWifa) ■ D S2 y j w j (t 2 ) ■ t 2 )0(s 2 , ti) ds 2 dsi dt 2 dti 

'U JU 

D S2 y j w i (t 1 )(f)(s2,t 1 )ds2dt 1 - / D si y i w j (t 2 )4>(s u t 2 )ds 1 dt 2: 



U Jo 



ti JO 



we have by Lemma 15 that 

E(5(^Wil[ tl , tl+l] )5(3 ; i^lfe,t J+1 ] 

= E / <l+1 / <J+1 y^j- • ^(sH(t) • 0(s,t) rftrfs + 0(A 4H+2 ), 

which shows the assertion. 



□ 



Now we finally determine the strong asymptotic behaviour of 



E 



JO 



\{B t -B\ n) )dt 



For similar calculations in the case that the weight function y is deterministic 
and the process B is stationary and behaves locally like a fractional Brownian 
motion, see e.g. [34] and [2]. 
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Proof of Proposition 8. By Lemma 16 we have 

2\ 1/2 



E 



y t (B t - B[ n) ) dt 



(n—ln—1 ft +i rt +i \ 

i=0 j=0 Jti Jt i J 



Thus, it remains to study the behaviour of 



n— 1 n— 1 



EE/ / Ey i y j -w i (s) Wj (t)- ( f)(s,t)dtds. 



Note that 

'•ij+l ftj+l 



U Jt 



Wi{S)Wj 



(t) ■ (f)(s, t) dtds = E f tl+1 Wi{s) dB s Wj(t) dB t 



by (11) in Section 2.1, and recall that 



fU+i ft i+1 i . . 

I Ms) dB s = J - (B u + B ti+1 ) - B s ds. 



Define 



e id ( Sl , s 2 ) = -E(B U + B ti+1 - 2B Sl )(B tj + B tj+1 - 2B S2 ) 
for Si G [ti, U+i], s 2 G [tj, tj+i], i,j — 0, . . . , n — 1. Thus we can write 

n-ln-l rt'+i /■*'+! 

EE// / 3 • ^(s)^(t) • 0(s,t) cftote 



i=0 j=0 ,y * 



n— 1 n— 1 



EE / / E ^ ■e hJ ( Su s 2 )ds 2 ds 1 . 

i=0 j=Q Jt i Jt - 1 



By straightforward calculations we obtain 

$ij( s l) s 2) = _ - (\U — tj\ 2H + \ti — tj + i\ 2H + \t i+ i — tj\ 2H + \t i+ i — tj + i\ 2H ^j 



+ 4 {\u - s 2 \ zn + \u+i - s 2 r + \tj - sir + \t j+ i - si 



2H 
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which simplifies in the case % = j to 

1 

4 
1 

4 
1 



9j,j(si, s 2 ) - —-\t j+1 - tj\ 2 1 



+ \ [\tj ~ S2\ 2H + \t j+ l - S 2 \ 2H + \tj - Sl \ 2H + \t j+1 ~ Si| 2 " 



2H 



(i) We first show that asymptotically the contribution of the off-diagonal terms 
to the error is negligible, i.e., 



/ / Vy t y r 6 hJ (s 1 ,s 2 )ds 2 ds 1 = o(n- 2H - 1 ). (A.2) 



|i-j|>log(n) 

Note that by symmetry 



22 Ey i y j -9 i j(s 1 ,s 2 )ds 2 ds 1 

,1>log(n) Jt * Jt J 

rU+i rt j+1 

= 2- / / vy i yj-diAs 1 ,s 2 )ds 2 d Sl 

. , / \ J to J tj 



j-j>log(n) 

To show (A.2) we will use fourth order Taylor expansions of suitable parts of 
9i,j{si, s 2 ). For this, the following will be very helpful: 

Let a > and e x , e y G { — 1, +1}- Then for x, y G [0, 1] such that a + e x x + e y y > 
0, the function 

f(x,y) = (a + e x x + e y y) 2H 
is well defined and we have 

°" f —(x,y) = K n -e k x e^ k -(a + e x x + e y y) 2H - n (A.3) 



{dxf{dy) 
with 

K n = 2H ■ (2H - 1) • . . . • (2H - n + 1). 
In what follows, set 

Ti,j = \U+i/ 2 — tj+l/2 1 = A|i- j\, 

and recall that 

Wj(Si) = Si - ti+l/2, Wj(s2) = S 2 - *j+l/2- 
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The first part of 6 i: j we study is 

/)(!)_ 1 fu + \2H , |. . \2H , i. + |2i? , U f |2tf\ 

"ij — ~~ g \J r i — r jl + l r « ~~ + ~~ l j\ "+" R+l ~~ J 

= -g (2r^f + |r u - A/2 - A/2| 2 " + |r w + A/2 + A/2| 2 ") . 

Since i > j, we obtain by applying (A. 3) with a = Tij, x = A/2 and y = A/2 
the expansion 

M) 1 2H K 2 2H-2 A 2 , „ 

with 

|p M | <c-|t,-t, +1 | 2H - 4 -A 4 . 
Since i — j — 1 > log(n), we have 

p itj = O (log(n) 4 - 2 "n- 4 ) . 



Now consider the second part of given by 

0?M, S2) = \ ((U - S 2 f H + (t i+1 - S 2 f H ) 

= \ ((ry - ^2) " A/2) 2 " + (r M - ^(s 2 ) + A/2) 2 ") . 



4 
1 

4 

Here we obtain by applying (A. 3) with (s 2 ) and y = A/2 



*0 = \- rt» + | • r 2 f - ■ ((* 2 - t J+1/2 ) 2 + A 2 /4) 



T i,f 1 ' ( S 2 - ^+1/2) 



2 

12 

^3 r 2i/-:i 

16 ' ' 



r i\j ' ( s 2 - ^+1/2) • A 
+ O (log(n) 4 - 2 "n- 4 ) . 

Note that clearly 

/ s 2 - tj+i/2 ds 2 = 0, / (s 2 - tj+i/2) ds 2 = 0, 

3 3 

so the third, fourth and fifth term of the above expansion vanish after inte- 
gration over [tj,tj +1 ] x [tj,tj + i\. For 

0ZHsi,S2) = \{(si-t i ) 2H + (8l-t S+1 )™) 

= \ (fa + Wi ( Sl ) + A/2) 2 " + ( TiJ + Wi( Sl ) - A/2) 2 ") 
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we obtain analogously 



1 2H , K 2 2H-2 
— • T- ■ + ' T- • 

2 *j 4 *j 
+ (log(n) 4 - 2i/ n- 4 ) , 



si - U +1/2 f + A 2 /4) + i/ij 



where u it j denotes the terms of the Taylor expansion, which contain odd powers 
of si— U+i/2 resp. S2—tj+i/2 and vanish after integration over [ti, U+i\ x [tj, tj+i\. 

Finally, we have to study 



Oij(s U S 2 ) = -^|S1 - S 2 | 2W = ~\Tij + Wi(si) - Wj(s2, 



\2H 



Here we obtain 



0Q(si, s 2 ) = -\ ■ r 2 f - *j ■ r 2 f - 2 • - U +1/2 f + (s 2 - t j+1/2 ) 2 ) + vy 

+ o (\og( n y 



where denotes again the terms, which vanish after integration. 
Summing up the above expansions for the parts of 9 it j yields 

Msi, s 2) = Vij + O (\og(n) 4 - 2H n- 4 ) . 
Therefore we obtain after integrating over [tj,t i+ i] x [tj,tj +1 ] the estimate 



■i+1 ftj+1 



U Jt 



6i,j(si,s 2 ) ds 2 dsi 



<c(log(n)) 4 -^.-. 



Since y is bounded, it finally follows 

ftj+1 



rn+i ftj+i 

E / / E^^- • 0ij(si, s 2 ) rfs 2 dsi 



|i-j'|>log(n) ' 

and hence we have shown (A. 2). 



< c- (log(n)) 



4-2i/ 



r?; J 



(ii) Now it remains to consider the summands with \i — j\ < log(n), i.e. the 
diagonal resp. near diagonal parts. For this, we need to compute the integrals 
over [ti,t i+ i] x [tj,tj + i] of the four parts of 9 it j. Note that by symmetry 



E 

|»-j'|<log(n) 



U+i rtj+i 

tj, J tj 



2" E 

0<i— j<log(n) 



E^j^j • Oi,j{si, s 2 ) ds 2 dsi 

U+l ftj+1 



ti J t 



E^j^j • Oi,j(si, s 2 ) ds 2 dsi 



ra-1 



+ E 



tj+1 ftj+1 



• 0jj( s l) S 2) ^S 2 dSl- 
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U+l ftj + 1 



For i > j we obtain 

Qij(si, s 2) ds 2 dsi 

— — ^- (\U — tj\ 2H + \U — tj + i\ 2H + — tj\ 2H + \t i+ i — tj + i\ 

= 4 • (21* " ^ + |< + 1 - J| M + |< - 1 - ;T) • (A.4) 
Moreover, we have 

rU+i ftj+i 



fH+l (2) 

/ / 9-j(s 1 ,s 2 )ds 2 ds 1 
= 4(2//+ 1) ~ tj ' 2H+1 + ~ ^I 2i?+1 ) 



4(2//+ 1) 

1 (li+i-^-n-i-ii-^). ' (a.5) 



4(2J3"+1)V 1 1 1 ' n 

rU+i ftj+i 



fH+l ftj + 1 /QN 

/ / Oij(si,s 2 )ds 2 dsi 

J £ j j 

= 4(2// + 1) ~~ ^'| 2H+1 + ~~ 



2H+1 



4(2//+ 1) 

i _ ir+ . r«) i (a.6) 



4(2//+ 1) v Jl 1 Jl / n 

and 

*»+i ftj+i _( 4 ) 



8ij(si,S2)ds2 dsi 

I; /+ / \2H+2 , I . . |2fT+2\ 

2(2//+l)(2// + 2) U^-^l ) 



2(2//+ 1) (2// + 2) 



— ^'+i| 2 ' f/+2 + \U+i — ^'| 2H+2 ) 



1 ( 2 |i - j| 2 ^+ 2 



2(2//+ 1) (2// + 2) 



N-i-jr +2 -N+i-jr +2 )--^. (a.7) 



Moreover, for % = j we have 



Oj,j(s!, s 2 ) ds 2 dsi 



u -it 



\2H + 1 (2// + 2)(2//+ 1) A) n 2H + 2 ' 



(A.8) 
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Hence, combining (A. 4), (A. 5), (A. 6) and (A. 7) we obtain for % > j that 

rU+i ftj+i 1 
/ / Oi,j(s 1 ,s 2 )ds 2 ds 1 = ■ ld(i,j) 

with 

MiJ) = ~ ■ (2|i -j\ 2H + l> + 1 ~0\ 2H + \i - 1 - il 2 ") 
In the case i = j we have by (A. 8) that 

rU+i rtj+i I 

J u It, 6i ^ Su S2) ds2 dSl = 2 ' 

with 

r = ^ 1 1 

2 2H+1 (2H + 2)(2H+l) 4' 

Hence it follows 

rU+i ft j+i 

E / / e Wj ' s s) rfs 2 dsi 



2H+2 



|»-j|<log(n) 



I ?^_..o 1 



0<j-i<log(n) ,l j=0 



Since by Lemma 13 it holds 

1 



E|^-^r<clogH 2 — 



for \k — l\ < log(ra), we have 

E f !+1 E W" • M s i> s 2) ds 2 ds! (A.9) 



|»-j|<log(n) 



0<i-i<log(n) j=0 
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Rearranging the terms on the right hand side of the above equation yields 

2. £ e^mmi^+ee^^ 1 



rl 2H+2 ' ^ 3 z rl 2H+2 
0<j-j<log(n) 11 j=0 11 

n— 1 j 



•C,-(Llog(n)J), 

i=o ™ 



where 



C,(r) = /C 2 + 2 £ /Ci(i,j) 

i=j+l 

for j = 0, . . . , n - 1, r e N. Note that 

2 £ /C 1 (i,j) = ~E2fc w + (fc+l) M + (fc-l) 2fl 
i=j+l 4 fc=l 



+ (2iJ + 1) 1 (2iJ + 2) £ ^ - (* + i) 2H+2 - {k - i) 2H+2 

= -W+i) 2H -+-r 2H +]-f:k 2H 

4 4 4 ti 

+ ^— i ((r + l) 2 ^ + r 2H ^-l) 

+ (r 2H+2 - (r + l) 2H+2 + 1 

(2#+l)(2# + 2) V K ' 

for r G N. Thus we have 

C,-(r)=Co(r) 

for j = 0, . . . , n — 1 and 



1, 1 ,, 

r — 

fe=l 



Cb(r) = ~(r + 1)^ + ^-5:^ 



+ ^ T r((-+ir +1 +^ +1 ) 

(2# + l)(2# + 2) V 1 ; / 



In the next step we show that 

lim Co(r) = -C(-2H). (A.10) 
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For this, we again apply Taylor expansions of suitable parts of C (r). We have 
1 ( r + 1) 2 "+ 2 = 1 r 2H+2 + — l —r 2H+l 



(2#+ 1)(2F + 2) V ' (2H+l)(2H + 2) 2H + 1 

+ \r 2H + f r™- 1 + o(l), 

' ( r + 1)2H+1 = _J_ r 2H+l + r 2if + ^ r 2H-l + o(1)) 



2H+V ' 2H+1 

hr + l) 2H = 1 -r 2H + ^r 2H - 1 + o(l). 
4 4 2 



Hence it follows 



Coir) = ~tf H + ^Vr 2 * +1 + + fr 2 ^ + 



Since 



see e.g. [1], we obtain (A. 10). 
This yields 

Jim n 1+2H J2 f l+1 • s 2 ) ds 2 ds x 



|i-j'|<log(n) 



n-1 



= jim Co(LiogWj) • E Eyj • (t s+1 - 1 3 ) 

3=0 

n-1 

= \C(-2H)\-Vm'ElSyZ-(t j+1 -t j ) 

3=0 

= \C(-2H)\. fEfldt, 
Jo 

which together with (A. 2) shows finally Proposition 8. 



□ 
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